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We calculate the B — » D^*> form factors in the heavy-quark and large-recoil limits in 
the perturbative QCD framework based on kr factorization theorem, assuming the hierachy 
Mb 3> M D {t) ^> A, with the B meson mass Mb, the meson mass M D (»), and the heavy 
meson and heavy quark mass difference A. The qualitative behavior of the light-cone D^*' 
meson wave function and the associated Sudakov resummation are derived. The leading-power 
contributions to the B — > form factors, characterized by the scale Ay/ Mb /M D (») , respect 
the heavy-quark symmetry. The next-to-leading-power corrections in 1/Mb and 1/M D (,) , char- 
acterized by a scale larger than are estimated to be less than 20%. The D^*> meson 
wave function is determined from the fit to the observed B —* D^*'lv decay spectrum, which 
can be employed to make predictions for nonleptonic decays, such as B — > D^*\(p). 



* E-mail: 
f E-mail: 

* E-mail: 



krmt@sci.toyama-u.ac.jp 

hnli@phys.sinica.edu.tw 

sanda@eken.phys.nagoya-u.ac.jp 



I. INTRODUCTION 



Recently, we have made theoretical progress in the perturbative QCD (PQCD) approach to heavy-to-light 
decays, such as the proof of kr factorization theorem [1], the construction of power counting rules for various 
topologies of decay amplitudes [2] , the derivation of kr and threshold resummations [3-5] , and the application to 
heavy baryon decays [6] and to three-body nonlcptonic decays [7] . Important dynamics in the decays B — > Kit 
and B — > tttt has been explored, including penguin enhancement and large CP asymmetries [2,8-10]. It is 
worthwhile to investigate to what extent this formalism can be generalized to charmful decays, such as the 
semileptonic decay B — > D^lv and the nonlcptonic decay B — > D^n(p). kr factorization theorem for the 
B — > form factors in the large-recoil region of the L>W meson can be proved following the procedure in 
[1], which are expressed as the convolution of hard amplitudes with B and £)W meson wave functions. A hard 
amplitude, being infrared-finite, is calculable in perturbation theory. The B and meson wave functions, 
collecting the infrared divergences in the decays, are not calculable but universal. 

The B — > transitions arc more complicated than the B — > tt ones, because they involve three scales: 
the B meson mass M B , the £>W meson mass M D (,) , and the heavy meson and heavy quark mass difference 
A = Mb — nib ~ M fl( «) — m c , where mb (m c ) is the 6 (c) quark mass, and A of order of the QCD scale Aqcd- 
There are several interesting topics in the large-recoil region of the B — > £)W transitions. (1) How to construct 
reasonable power counting rules for these decays with the three scales? (2) What is the qualitative behavior 
of the wave function for an energetic meson? Is it dominated by soft dynamics the same as a B meson 
wave function, or by collinear dynamics the same as a pion wave function? (3) How different are the end-point 
singularities in the B — ► Z)W form factors from those in the B — > 7r ones? (4) what is the effect of Sudakov 
(kx and threshold) resummation in the B — > form factors? (5) How are the relations for the 2? — > Z)(*) 
form factors in the heavy-quark limit modified by the finite B and meson masses? (6) What is the hard 
scale for the B — ► £)(*) transitions? Is it large enough to make sense out of perturbative evaluation of hard 
amplitudes? These questions will be answered in this work. 

We attempt to develop PQCD formalism for the B — ► D^*Hv decays in the heavy-quark and large-recoil limits 
based on kx factorization theorem [11,12]. We argue that the following hierachy must be postulated: 

M B > M DW > A . (1) 

The relation Mb *> M D (.) justifies the perturbative analysis of the B — > £)(*) form factors at large recoil and 
the definition of light-cone 

£>(*) meson wave functions. The relation M D (,) 3> A justifies the power expansion in 
the parameter A/M D (,>. We shall calculate the B — > form factors as double expansions in M D (.) /Ms and 
in A/M D (,). The small ratio A/Mb = (M D (,)/Mb)(A/M D (»)) is regarded as being of higher power. Whether 
the hierachy is reasonable can be examined by the convergence of high-order corrections in the strong coupling 
constant a s and of higher-power corrections in I/Mb and in 1/M D {»). Note that Eq. (1) differs from the 
small- velocity limit considered in [13], where the ratio M D (*) /Mb is treated as being of 0(1). 

Under the hierachy, the wave function for an energetic Z)W meson absorbs collinear dynamics, but with 
the c quark line being eikonalized. That is, its definition is a mixture of those for a B meson dominated by 
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soft dynamics and for a pion dominated by collinear dynamics. We examine the behavior of the heavy meson 
wave functions in the heavy-quark and large-recoil limits. For A/Mb, A/M fl (.) <C 1, only a single B meson 
wave function 4>b{x) and a single Z}W meson wave function <fi D (,)(x) arc involved in the B — > Z)W form 
factors, x being the momentum fraction associated with the light spectator quark. Equations of motion for 
the relevant nonlocal matrix elements imply that 4>b{x) and <p D (,){x) exhibit maxima at x ~ A/ Mb and at 
x ~ A/M D (»), respectively. Since PQCD is reliable in the large-recoil region, it is powerful for the study of 
two-body nonlcptonic B meson decays. After determining the £)W meson wave function (the B meson and 
pion wave functions have beed fixed in our previous works), we can make predictions for the B — > D^tt(p) 
decays. Preliminary results for the B — ► Dn modes will be presented in Sec. IV. 

Similar to the heavy-to-light transitions, the B — > form factors also suffer singularities from the end- 
point region with a momentum fraction x — > in collinear factorization theorem. When the end-point region is 
important, the parton transverse momenta kr are not negligible, and kx factorization is a tool more appropriate 
than collinear factorization. It has been proved that predictions for a physical quantity derived from kr factor- 
ization are gauge-invariant [1] . Because of the inclusion of parton transverse degrees of freedom, the large double 
logarithmic corrections a s In 2 kr appear and should be summed to all orders. It turns out that the resultant 
Sudakov factor for an energetic £)M meson is similar to that for a B meson. Including the Sudakov effects from 
kr resummation and from threshold resummation for hard amplitudes [5,14], the end-point singularities do not 
exist, and soft contributions can be suppressed effectively. 

We shall identify the leading-power and next-to-leading-power (down by 1 /Mb or by 1/M D (,)) contributions 
to the B — > £)M form factors, which are equivalent to the expansion in heavy quark effective theory [15,16]. 
It will be shown that the leading-power factorization formulas satisfy the heavy-quark relations defined in 
terms of the Isgur-Wise (IW) function [15]. These contributions are characterized by the scale A^/Mb/M D (,), 
indicating that the applicability of PQCD to the B — > form factors may be marginal. The next-to-leading- 
power corrections to the heavy-quark relations, characterized by a scale larger than \/AMb, can be evaluated 
more reliablly. These corrections amount only up to 20% of the leading contribution, indicating that the power 
expansion in 1/Mb and in 1/M D («) works well. It should be stressed that the above percentage is only indicative, 
since we have not yet been able to explore the complete next-to-leading-power sources with the current poor 
knowledge of nonperturbative inputs. 

In Sec. II we define kinematics and explain kr factorization theorem for the B — ► D^*Hv decay. The power 
counting rules are constructed. In Sec. Ill we discuss the behavior of the B and meson wave functions, and 
perform kr resummation associated with an energetic meson. The leading-power and next-to-leading-power 
contributions to the B — ► £)W form factors are calculated in Sec. IV. Section V is the conclusion. 

II. KINEMATICS AND FACTORIZATION 

We discuss kinematics of the B — > D^*Hf decay in the large- recoil region. The B meson momentum Pi and 
the meson momentum P 2 are chosen, in the light-cone coordinates, as 
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Mr A^Mr 

tt = -^(l,l,Or), P2 = ^^(v + ,V-,0t), (2) 



with the ratio A*^ = M D (») /Mb- The factors r/ ± — r\ ± i s defined in terms of the velocity transfer 

r\ = v\ ■ V2 with v\ — Pi /Mb and t>2 = P 2 /M D (*). The longitudinal polarization vector cl and the transverse 
polarization vectors ct of the D* meson are then given by 

e L = -L(n+,-ri-,0T) , er = (0,0,1). (3) 

The partons involved in hadron wave functions are close to the mass shell [17]. Assume that the heavy meson, 
the heavy quark, and the light spectator quark carry the momenta Ph, Pq = Ph — k, and k, respectively. The 
on-shell conditions lead to 

k 2 = 0(A 2 ) , (4) 
P 2 - m% = M 2 H - m% -2P H -k = 0(A 2 ) , (5) 

Mh (juq) being the heavy meson (quark) mass. For the B meson, we have 2P\ ■ k\ — 2Msk\ ~ M^ — mf } = 
2Mb A from Eq. (5), namely, kf ~ A. The order of magnitude of the light spectator momentum is then, from 
Eq. (4), 

K ~ (A, A, A) . (6) 
For the £)(*) meson, we have 2P 2 ■ k 2 = 2P 2 k 2 + 2P 2 k 2 ~ 2M D( *)A, which implies 

2 \M DM M B J y ' 

for r] + ~ I/A** 1 and r\~ <~ r^. With the above parton momenta, the exchanged gluon in the lowest-order 
diagrams is off-shell by 

(*i - fe) 2 - -T^A 2 , (8) 

M D (.) 

which is identified as the characteristic scale of the hard amplitudes. To have a meaningful PQCD formalism, 
the large-recoil limit in Eq. (1), Mb ^> M D (,), is necessary. 

We have proved fc-r factorization theorem for the semileptonic decay B — ► 7rZf [1]. Soft divergences from the 
region of a loop momentum Z, where all its components are of 0(A), are absorbed into a light-cone B meson 
wave function. Collinear divergences from the region with I parallel to the pion momentum in the plus direction, 
whose components scale like 

l"~ (M B , A 2 /M b , A), (9) 

are absorbed into a pion wave function. The above meson wave functions, defined as nonlocal hadronic matrix 
elements, are gauge-invariant and universal. 

fcx factorization theorem for the semileptonic decay B — > D^*Hv is similar. In the limit Mb 3> M D (,) we 
have k 2 ^> k 2 T 3> k 2 from Eq. (7), indicating that the £)W meson wave function is dominated by collinear 
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dynamics. The leading infrared divergences in this decay are then classfied as being soft, if a loop momentum I 
vanishes like <~ (A, A, A), and as being collinear, if I scales like k 2 in Eq. (7). The former (latter) are collected 
into a light-cone B (P/*-*) meson wave function. The collinear gluons defined by Eq. (9) do not lead to infrared 
divergences. 

Though the P/*) meson wave function absorbs the collinear configuration, similar to the pion wave function, 
the heavy-quark expansion applies to the c quark in the same way as to the b quark in a B meson. This is the 
reason we claim that the energetic meson dynamics is a mixture of those of the B meson at rest and of the 
energetic pion. Because P 2 • Z ~ M D (,)A is much larger than I 2 ~ A 2 according to Eq. (1), we have the eikonal 
approximation, 

Pi- fa+ /- 



n a c(P 2 - k 2 ) 



F c(P 2 -fc 2 ) 



(10) 



(P 2 - k 2 + I) 2 - m 2 c ' ' v 2 -V 

where c(P 2 — k 2 ) is the c quark spinor, and the factor v 2 j(v 2 ■ I) the Feynman rule for a rescaled c quark field. 
The physics involved in the above approximation is that the kinematics of the spectator quark and of the c 
quark is dramatically different in the limit M D (») — > oo. Hence, a gluon moving parallel to k 2 can not resolve 
the details of the c quark, from which its dynamics decouples. 

According to kr factorization theorem [1], the light spectator momenta k\ in the B meson and k 2 in the £)W 
meson are parametrized as 



in Mb i, 

M ( 0, Zi—/=-, kiT 



k 2 



where the momentum fractions x\ and x 2 have the orders of magnitude, 



x\ ~ A/Mb , x 2 ~ A/M D 



(.) 



(11) 



(12) 



The smallest component k 2 has been dropped. The neglect of fc+ is due to its absence in the hard amplitudes 
shown below. 





q q 

FIG. 1. Lowest-order diagrams for the B — > form factors. 



The lowest-order diagrams for the B — > form factors are displayed in Fig. 1. The factorization formula 
is written as 



(D^(P 2 )\bT fl c\B(P 1 )) =g 2 C F N c I d Xl dx 2 d 2 k 1± d: 2 k 2± 



I- 



v2,, .2, dz+d 2 z ± dy d 2 y ± 



(2tt) 3 (2tt) 3 

xe- ife2 ^(£)W(P 2 )|J 7 ( 2/ )c /3 (O)|O)e ifel -(O|6 Q (O)^(0)|v3(P 1 )) Hf- aS 



(13) 



with the hard amplitude, 
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b«V 5 - 



1 



(fc 2 



h) 2 
1 



■ /^2~ /Pi + rah 

(-Pi - k 2 y - 



2^ 



" l 6 



a/3 



a/3 



(fe-^) 2 r M (^-*i) 2 -"»2 r j ' (14) 

for T M = 7^ or 7^75. It is obvious that the large component picks only up the component fcf in the 
denominators of the internal particle propagators. The first and second terms in behave like M^ (t) /(A 3 Af|) 
and M D (,) /(A 3 Mb), respectively. Therefore, for a leading-power formalism, we keep 0(1) coefficients of the 
first term, and 0(A*^) coefficients (0(1) coefficients are absent) of the second term. 



III. HEAVY MESON WAVE FUNCTIONS 



In this section we discuss the qualitative behavior of the B, D and D* meson wave functions in the heavy-quark 
and large-recoil limits, and derive kr resummation associated with the meson. 



A. B Meson Wave Functions 



According to [1,18,19], the two leading-twist B meson wave functions defined via the nonlocal matrix element, 



/ 



d w 
(2^ 



e ik - w (0\b a (0)d s (w)\B°(P)) = - 



(/P + M B ) 75 



^U+(fc) + 4<Mfc) 



V2 



a 5 



/2N, 



with the dimensionless vectors n = (1, 0, Ot) and n = (0, 1, Ot) on the light cone, and the wave functions, 



J B ' 



(15) 
(16) 

(17) 



We have shown that the contribution from <j) B starts from the next-to-leading-power A/M B . This contribu- 
tion, which may be numerically relevant [20], should be included together with other next-to-leading-power 
contributions in order to form a complete analysis. On this point, our opinion is contrary to that in [21]. 

The investigation based on equations of motion [22] shows that the distribution amplitude (j>^(x) — 
J d 2 kT(j) B (x,kT) vanishes at the end points of the momentum fraction x = k~ / P~ — > 0, 1. Hence, we adopt 
the model in the impact parameter b space [8], 



b B (x, b) — N B x 2 {\ — x) 2 exp 



1 ( xM B 

2 I co B 



oo 2 B b 2 



(18) 



where the shape parameter has been determined as ujb — 0.4 GeV. The normalization constant Ng is related 
to the decay constant fs through 

Sb 



I 



dx(t> B {x, 6 = 0) 



2J2N r 



(19) 



It is easy to find that 4>b in Eq- (18) has a maximum at x ~ A/Mb as claimed in Sec. II. 
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B. D Meson Wave Functions 



Consider the nonlocal matrix elements associated with the D meson, 



(0|c(y) 7 57^HI^"(^)> = -ifnP,* f dxe^-^-^ cjf D {x) 

Jo 

--JuMl^- [ 1 dxe-« p - k >y- ik -*>g D (x) , (20) 
1 " ' z Jo 

(0\c(y) l5 d(w)\D-(P)) = -if D m f foe-HP-V-v-*-* ^ D ( x ) , (21) 

Jo 



M 2 

(0\c(y) l5 a^ u d(w)\D-(P)) = --f D m ( 1 §■ ) (P„*„ - P vZfi ) 

1Tln 



x f dxe-^ p - k > y -* k - w (j) a D (x) , 
Jo 



(22) 



with z = y — w and the D meson decay constant fo ■ The light spectator d quark carries the momentum k with 
the momentum fraction x = k + / P + and the c quark carries the momentum P — k. In the heavy-quark limit 
we have 

M 2 

m = 2_ = Md + 0(A) , (23) 

m c + m d 

implying that the contribution from the distribution amplitude (\f D is suppressed by 0(A/Md) compared to 
those from <p v D and <jf D . The distribution amplitude go, appearing at 0(r 2 ), is negligible. 
Rewrite the pseudo-tensor matrix element as 

(0\c(y) l5 a^ v d(w)\D-(P)) = i(0\c(y) l5 j^d(w)\D- (P)) - ig^(0\c(y) l5 d(w)\D- (P)) , (24) 

and differentiate both sides with respect to w and y. The differentiation on the left-hand side gives a result 
suppressed by 0(A/M D ). The relations 

J dxx^ p D (x)e- l(p - k >y- lk - w = 0(A/M D ) , (25) 

J dxW D {x) - r D (x)]e-^ p - k >y- ik - w - 0(A/M D ) , (26) 

arise from the differentiation with respect to w v for fi = — and to y^ for v = + , respectively. Equation (25) 
states that the distribution amplitude 4> P D possesses a maximum at x ~ A/Mn- Equation (26) states that the 
moments of <j? D and <p v D differ by 0(A/Mr>) (they have the same normalizations). 

Neglecting the 0{A/Md) difference according to Eq. (1), only a single D meson wave function is involved in 
the evaluation of the B — > D form factors, 

/ ^e ik ' W (QM^M\D-(P)) = --^[(JP + MoYf^Dix) , (27) 

where the distribution amplitude, 

* D - Sr/° = Ar/° ■ < 28 > 
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satisfies the normalization, 



fn 



/ dxd> D (x) = 
Jo VDK ' 2^2N C 

For the purpose of numerical estimate, we adopt the simple model, 



(29) 



(30) 



The free shape parameter Cd is expected to take a value, such that (fro has a maximum at x ~ A/M D ~ 0.3. 
We do not consider the intrinsic 6 dependence of the D meson wave function, which can be introduced along 
with more free parameters. Note that Eq. (30) differs from the one of the Gaussian form proposed in [23]. 



C. D* Meson Wave Functions 



The information of the D* meson distribution amplitudes is extracted from equations of motions for the 
nonlocal matrix elements, 

t-i 



\c(y)^d(w)\D*-(P,e)) = f D .M D 



P„ 



dxe-^-^y-^^ix) 



+e Tfl / dxe-^ p - k > y - ik - w g { l\x) 
Jo 

jf dxe^yy-^ 93 (x) 



(31) 



(Q\c(y)a^d(w)\D*-(P,e))=if^ 



JO 



dxe 



i(P—k)-y — ik-w 



(j>j_(x) 



+(P ll z„ - PvZ^j^^Ml, [ dxe-* p - k >y- ik - w hf{x) 
(/ ' z ) Jo 

1 M 2 f 1 

+ ( eT ^-e Tv z^— ^ / dxe-^ p - k >y- ik - w h 3 (x) , (32) 
z " • z Jo 

(0\c(y)d(w)\D*-(P,e)) = -\ (/£. f D * ^j^-) e • zM D * £ dxe^-^-^h^x) , (33) 

(0\c(yh 5 ^d(w)\D*-(P,e)) = -\ (j D * fl* ™^™* ) M D .e™Pe Tv P a z p J* dxe-^-^-^g^x) , (34) 

where the D* meson decay constant Jd* (/^, ) is associated with the longitudinal (transverse) polarization. 
In the heavy-quark limit we have 



Id* - Id 



m c + m d 



Jd* - fo* 



T m c + m d 



0(A/M D *) . 



Md* " JU Md* 
Hence, the contributions from the various distribution amplitudes are characterized by the powers, 

0,1 (x), 4>±(x) : 0(1), 
<#>(*), hf{x): 0(r*), 
g 3 (x), h 3 (x): 0(r* 2 ), 



(35) 



^ s) (x), g™(x): 0(A/M B ) 



» 



(36) 
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To the current accuracy, we shall consider the distribution amplitudes <j)\\ (x) and hfjp (x) for the longitudinal 
polarization, and (j)j_(x) and (x) for the transverse polarization of the D* meson. 
Rewrite the tensor matrix element as 

(0\c(y)a^d(w)\D*-(P)) = i{0\c(y)^ v d(w)\D*- (P)) - ig^{0\c(y)d(w)\D*- (P)) , (37) 

and differentiate both sides with respect to w and y. The relations, 

J dxxh { ^{x)e- l{p - k) - y - lk - w = 0{A/M D .) , (38) 

J dxx[<t> ± (x) + h 3 (x)]e- t( - p ~ k >y- tk - w = 0(A/M D .) , (39) 

J dx[(f,\\(x) - hf{x)]e-« p - k yy-^ w = 0{A/M D ,) , (40) 

J dx[^{x) - g 3 (x)]e-^ p - k >y- ik - w = 0(A/M D ,) , (41) 

J dx[<t>j_(x) + h 3 (x) - 2g i l ) (x)}e-^ p - k >y- lk - w = 0(A/M D .) , (42) 

come from the derivatives with respect to w v for /i = — , to w v for /i =_L, to for v = — , to y^ for v = +, and 
to y^ for v =_L, respectively. Equations (38) and (39) indicate that the D* meson distribution amplitudes have 
maxima at x <~ A/Mo*. Equations (40) and (41) state that 0||, h^' and g 3 are identical up to corrections of 
0(A/Md«). Similarly, </>_l, h 3 and g^ are also identical up to corrections of 0{A/Md*) from Eq. (42). 
Neglecting the 0(A/Md*) difference, we consider the structure for a D* meson, 

/ ^^ W mm^w)\D*-{P)) = --^=[{P + M D ,) ^ L D ,(x) + (/P + M D ,) £t4>1*{x)U > (43) 
with the definitions, 

- i7w/« = Sm h> " ■ (44) 

* " Sm^ " J^ 5 -" ' (45) 

The £>* meson distribution amplitudes satisfy the normalizations, 



l r i f 



dx^.(x)= I dxcp 1 D .(x) = ^=, (46) 



/o Jo 

where we have assumed fr>* = fp. . Note that equations of motion do not relate <j>p* and 4>d» ■ In this work we 
shall simply adopt the same model, 

- <^,(x) = -JL^.^l - x)[l + C D ,(1 - 2x)\ . (47) 

Similarly, the free shape parameter Co* is expected to take a value, such that <po* has a maximum at x ~ 
A/Afo. - 0.3. 
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D. Sudakov Resuraraation 



Radiative corrections to the meson wave functions and to the hard amplitudes generate double logarithms 
from the overlap of collinear and soft enhancements. The double logarithmic corrections a s In 2 kr to the heavy 
and light meson wave functions and their Sudakov resummation have been analyzed in [3]. The property of 
the £>W meson wave function is special, since its dynamics is a mixture of the soft one in the B meson wave 
function and the collinear one in the pion wave function. In this section we derive kr resummation for the 
6-dependent D meson wave function, 



b D (x,b) 



I 



dw Jk+Z 

2tt 



(0|c„ 2 (0)^-Pcxp 



•L 



ig / ds ■ A(s) 



d(w)\D-(P 2 )) , 



(48) 



with the coordinate w = (0, w~ , b) . The path for the Wilson link is composed of three pieces: from to oo along 
the direction of n, from oo to oo + b, and from oo + b back to w along the direction of — n [1]. The derivation 
for the D* meson is the same. The meson distribution amplitude c/) D (,)(x) discussed in the previous two 
subsections is regarded as the initial condition of the Sudakov evolution in b. 



(a) 



D 



(b) (c) 
FIG. 2. Radiative corrections to the D meson wave function in the axial gauge. 



In the axial gauge the double logarithms appear in the two-particle reducible corrections, such as Figs. 2(b) 
and 2(c). Figure 2(a) gives only single soft logarithms. To implement the resummation technique, we allow the 
direction n of the Wilson line to vary away from the light cone. Because of the rescaled c quark field, P 2 does 
not lead to a large scale, and the only large scale is k 2 . Due to the scale invariance in n and in v 2 , k 2 appears 
in the ratios (k 2 ■ n) 2 /n 2 and (k 2 ■ v 2 ) 2 /v 2 . However, the second ratio is of 0(r 2 ) compared to the first one, and 
negligible. Therefore, 4>d depends only on the single large scale (k\ ■ n) 2 /n 2 . The rest part of the derivation 
then follows that in [3]. The Sudakov factor from kx resummation is given by 



exp[-S DM ((j,)] = exp 



-s(k+,b)-2 



Ji/b M 



(M)) 



(49) 



with the quark anomalous dimension 7 = —a s /-K. For the explicit expression of the Sudakov exponent s, refer 
to [8] . It is found that Eq. (49) has the same functional form as the Sudakov factor for the B meson. 
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The double logarithms a s In 2 x produced by the radiative corrections to the hard amplitudes are the same 
as in the B — > ir,p decays at leading power in 1/Mb and in 1/M D <„). [14]. Threshold resummation of these 
logarithms leads to 

W' W ^ 1 '-' 1 ' (50) 

with the constant c = 0.3 ~ 0.4. The factor St(x 2 ) (St(xi)), associated with the first (second) term of in 
Eq. (14), suppresses the end-point region with x 2 — > {x\ — > 0). In the numerical study below we shall adopt 
c= 0.35. 



IV. B-*D,D* FORM FACTORS 

The B — > transitions are defined by the matrix elements, 



<£>(P 2 )|5(0) 7 „c(0)|B(Pi)) = V^K+WK + ^U^WK-^] , 



(D*(P 2 ,e*)|6(0) 7 ^75c(0)|£(Pi)) = ^M B M D , [Ui{v){v + 1)< - • - £«(»7)e* ■ «i«2 M ] 



( J D*(P 2 ,e*)|K0) 7M c(0)| J B(P 1 )) = i y/M B M D ^ v elv 2a v lf} . (51) 
The form factors £+, ^j, £a 2 j £a 3) and £y satisfy the relations in the heavy-quark limit, 

i+ = £,v = U = U = £, £- = U = , (52) 

where £ is the IW function [15]. 

We write the form factors as the sum of the leading-power and next-to- leading-power contributions, 

6=& (0) +tf\ (53) 
for i — +, — , A\, A 2l A3, and V. The leading-power factorization formulas are given by 

q_ 0) = 16TTC F VrM B J dxidx 2 J b\db\b2<lb2(t>B{xi,bi)<j)D{x2) 

x E{t^)h{x 1 ,x 2 ,b 1 ,b 2 ) + rE{t^)h{x 2 ,x 1 ,b 2 M)\ , (54) 

^ A2 = , (55) 
€ai,A3,v = IGkCfVt^Mb J dxidx 2 J bidbib 2 db 2 (/)B(xi, h)4> D * (x 2 ) 

x E(t^)h(x 1 ,x 2 ,b 1 ,b 2 ) + r*E{t^)h{x 2 ,x 1 ,b 2 ,b 1 )] , (56) 

with the color factor Cf = 4/3. Obviously, the above expressions obey the heavy-quark relations in Eq. (52), 
if we assume 

M D = M D , , <j> D = cj) D , . (57) 

The hard function is written as 
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h(x 1 ,x 2 ,b 1 ,b 2 ) = K (y/xiXirWri+MBlh) S t (x 2 ) 

x Ofa - b 2 )K (y X2 r^r 1 +M B b^ I (v 'x 2 r(*)r,+M B b 3 ) 
+8(b 2 - h)K (y/x 2 r(*)ri+MBb2) h ( vW ( * VA'/ B &i) . 

In the evolution factor, 

E(t) = a,(t) cxp[-S* B (t) - S DM (t)} , 



(58) 



(59) 



we keep the Sudakov factor associated with the B meson, and allow the behavior of the B meson wave function 
to determine whether this effect is important. For the model <f>B(x,b) in Eq. (18), the Sudakov effect is not 
important because of x <~ A/Mb- The hard scales t are defined as 



= max {^x 2 A*)r]+M B , l/6i, l/b 2 ) 
= max (V'VVMii, 1/6!, I/62) 



(60) 



Equations (54) and (56) contain at most the logarithmic end-point singularities in the collinear factorization, 
which are weaker than the linear singularities in the B — > 7r form factors. We conclude that Sudakov resumma- 
tion is also crucial for the B — ► transitions. 
The next-to-leading-power corrections are given by 



j dxidx 2 J bidbib 2 db 2 4>B(xi,bi)<j) D (x 2 ) 



( V + - !)(„+ - 2) 



rx 2 E(t (1) )h(x 1 , x 2 , bi, b 2 ) + x 1 E(t {2) )h(x 2 , xi, b 2 , bi) 



e = 



j dx\dx 2 J bidbib 2 db 2 4>B{x\,bi)(j)D{x 2 ) 



(r?+ + l)(r7+ - 2) 



E{t {1) )h{x u x 2 , 61, 6 2 ) - x 1 E(t (2) )h(x 2 , x u 6 2 , 61) 



f (i) 



SA2 



J dx\dx 2 J bidbib 2 db 2 (f) B {x\,bi)(j) D *{x 2 ) 



rj+ - 2 



- r*x 2 E(t^)h(x u x 2 , b u b 2 ) ±-E(tW)h(x 2 , Xl ,b 2 , h) 

T)+l 7] + 1 

-16nCFVr*M B J dx x dx 2 J bxdbxb 2 db 2 ^ B {xx, b\)$r>* (x 2 ) 



rj^xi 



-.E{t^ ) )h(x 2 ,x 1 ,b 2 ,b 1 ) , 



J dx\dx 2 J bidbib 2 db 2 4>B{x\,bi)(j)D'{x 2 ,b 2 ) 



-.r*x 2 E{t (1 ^)h(x u x 2l b 1 ,b 2 ) 



Xi 



-.E(t (2 ^)h(x 2 ,x 1 ,b 2l b 1 ) 



(61) 



(62) 



(63) 



(64) 



(65) 



. yV - 1 V 7 / 2 1 

whose hard amplitudes are consistent with those obtained in [24-26]. The additional powers in x 2 and in x\ 
provide stronger suppression in the end-point region with x\, x 2 — ► 0. Hence, the characteristic hard scales of 



the corresponding terms increase to y KMb and to Mb \/ A/M D (,-, , respectively. 
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Consider the expansion of the currents to 0(l/rrib) and 0(l/m c ) [27], 

bT l c w fe^Tc^ + — !— b Vl Tii p 2 c V2 - -^—b Vl Tii pxc V2 H , (66) 

where • • • stand for the terms, which arc further suppressed by a s . Compared to the right-hand side of Eq. (66), 
£^ and the terms proportional to x 2 (x\) in are identified as the first term and the second (third) term. 
We do not distinguish Mb and mi,, and M B (,), m c and mo, and employ only one single wave function for the 
B and £>W 

mesons. The differences of the above quantities are also the sources of 1/m corrections. However, 
their estimation requires more information of nonperturbative inputs, and can not be performed in this work. 
Equations (61)-(65) will be employed to obtain an indication of the order of magnitude of next-to-lcading-power 
corrections to the B — > form factors. 




(a) (b) 
FIG. 3. (a) Contribution to at the maximal recoil from the different ranges of a s /iv. (b) Contributions to £ 
(black) and to £^ (gray) at the maximal recoil from the different ranges of a s /n. 
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0.25 
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0.15 
0.1 

0.05 



1.35 
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1 . 45 



1.5 



v . v' 



(b) 



FIG. 4. (a) Dependence of (square) and of ^V^^ (triangle) on the velocity transfer, (b) Dependence of 

-^Ai/^Ai (square), of -£%1/£a1 (triangle), and of Ca^v/^Ai ( star ) on the velocity transfer. 

It is observed from Fig. 3(a) that most of the contribution to the form factor comes from the range of 
a s /-K < 0.3, implying that the applicability of PQCD to the B — > form factors is acceptable, and not worse 
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than that to the B — > ir ones [4]. This is attributed to the fact that the hard scales A-^/Mb/M^,) and \/KMb 
in the two cases do not differ very much. The applicability improves for the next-to-leading-power contributions 
as shown in Fig. 3(b): most of them arise from a s /Tr < 0.2. We emphasize that the above percentage analysis is 
only indicative, and that the convergence of higher-order corrections needs to be justified by explicit calculation. 
We estimate from Fig. 4 that next-to-leading-power contributions are less than 20% of the leading one. In fact, 
they are small except , implying that the power expansion in A/Mb and in A./M D (») is reliable. Our results 
are smaller than those obtained from QCD sum rules at large recoil [28]. 



form factor 



form factor 




- — • • 




(a) (b) 
FIG. 5. (a) [(b)] ^ as a function of the velocity transfer from the B — > D^*'lv decay. The solid lines represent the 
central values, the dashed (dot-dashed) lines give the bounds from the linear (quadratic) fits. The circles correspond to 
C D (,) — 0.5, 0.7, and 0.9 from bottom to top. 

The next step is to determine the meson distribution amplitudes, i.e., the free parameters C D (*), by 
fitting the leading-power PQCD predictions to the measured decay spectra at large recoil [29-31]. The IW 
function extracted from the B — ► D^Hv decay is parametrized as 



= F DW (1)[1 - p 2 DM ( V - 1) + c DM (v - I) 2 + 0(( V - l) 3 )] 



(67) 



with the factors F D (1) = 0.98 ±0.07 [32] and F D ,(1) = 0.913 ±0.042 [33]. The above values are consistent with 
those derived from lattice calculations [34]. The linear and quadratic fits give [29,30] 



Pn 
P% 



0.69 ±0.14, c D = 



0.81 ±0.12, cd« 



-1+0.42 
-0.15 ) 



00+ 59 
u - uu -0.00 ! 



(68) 



respectively. Choosing the decay constants fs = 190 MeV and fo = /u* = 240 MeV, we find that Co ~ 
Cd* = 0.7 leads to an excellent agreement with the data at large recoil as exhibited in Fig. 5. For these values, 
the corresponding Z}W meson distribution amplitude exhibits a maximum at i ~ 0.36, consistent with our 
expectation. The rough equality of Cd and Cd* hints that the heavy-quark symmetry holds well. 

Note that our aim is not to extract the Cabbibo-Kobayashi-Maskawa matrix element \V c b\ from experimental 
data. This extraction is best done in the zero-recoil region, where the heavy-quark symmetry defines unam- 
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biguously the normalization of the B — ► £)(*) transition form factors. It has been emphasized that the PQCD 
formalism is reliable at large recoil, and appropriate for two-body nonlcptonic decays. Therefore, one of the 
purposes of this work is to determine the unknown meson wave function. The B meson and pion wave 
functions have been fixed already in the literature. With these meson wave functions being available, we are 
able to predict the branching ratios of two-body nonleptonic charmful decays, such as B — > D^tt(p). Our 
predictions for the B — > Dir branching ratios [35], 



B(B~ - 




) - 5.5 x 1CT 3 


B(B° - 


■+ D+ir- 


) - 2.8 x 10" 3 


B(B° - 


- D\°] 


- 2.6 x 10~ 4 , 



(69) 

are in agreement with experimental data [36-38] . The above results correspond to the phenomenological coef- 
ficients oti and <22 [39] with the ratio |a2|/ai ~ 0.5 and the phase —57° of 0,2 relative to a\. The point is, from 
the view point of PQCD, that the phase is of short distance and generated from hard amplitudes. This is con- 
trary to the conclusion drawn from naive factorization [39-42]: the phase comes from long-distance final-state 
interaction. 

If the meson decay constant is known from, for example, lattice QCD calculation, it is then possible to 
extract the matrix element |V^i,| from the measured scmileptonic decay spectra at large recoil using the PQCD 
formalism. The experimental data of the product [|14(,|^(r/)] e xp for the B — > Dlv mode [29] are listed in Table I. 
The region with the large velocity transfer r\ > 1.35 is regarded as the one, where PQCD analyses are reliable. 
We compute the following %-square as a function of the two parameters, \V c b\ and the shape parameter Cd 
(iV para - 2): 

2 



x 2 = - y, 

-t*data ^para ^ 



|Kb|£te)]c XP HK 6 |d%fc^) 



07c 



(70) 



where A^data = 4, are assumed to be 10% of the data (considering only the systematic errors for illustration), 
and has been defined in Eq. (54). The band in Fig. 6 stands for the allowed |V c b|-Cj3 range with x 2 < 1. It 
is found that taking Cd = 0.7, V c b = (0.035 <~ 0.0036)(240 MeV / fo) is consistent with the value extracted in 
[29] from the zero-recoil data. Choosing V c b = 0.04, Cd = 0.4 <~ 0.5 is allowed. 
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0.03 0.035 0.04 0.045 0.05 



|Vcb| 



FIG. 6. The x 2 value as functions of \V c b\ and Cd with the shaded region corrsponding to \ 2 < 1. 



k 


1 


2 


3 


4 


V 


1.39 


1.45 


1.51 


1.57 


[\V cb \ari)U P 


0.028 


0.030 


0.024 


0.024 



TABLE I. Experimantal data of |V c fc|$(rj) for the B — > Dlv decay. 
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V. CONCLUSION 



In this paper we have developed the PQCD formalism for the B — ► transitions in the heavy-quark and 
large-recoil limits based on kr factorization theorem. The reasonable power counting rules for these decays 
with the three scales Mb, M D {,) and A have been constructed following the hierachy in Eq. (1). Under this 
hierachy, only a single B meson wave function and a single Z?W meson wave function arc involved, which 
possess maxima at the spectator momentum fractions x ~ A/Mb and x <~ A/M D( ,), respectively. Dynamics 
of an energetic meson is the mixture of those of a B meson at rest and of an energetic light meson: it 
absorbs collinear divergences but the heavy-quark expansion applies to the c quark. The Sudakov factor from 
kr resummation for an energetic meson is similar to that associated with a B meson. The end-point 
singularities, being logarithmic in collinear factorization theorem, do not exist in fcr factorization theorem. 
Including also the Sudakov effect from threshold resummation for hard amplitudes, the PQCD approach to the 
B — > transitions becomes more reliable. 

The factorization formulas for the B — > form factors have been expressed as the sum of leading-power 
and next-to-leading-power contributions, which is equivalent to the heavy-quark expansion in both 1/mt and 
l/m c . The leading-power formulas, respecting the heavy-quark symmetry, are identified as the IW function. 
This contribution, characterized by the scale K-J Mb/M D {») , is calculable marginally in PQCD. The next-to- 
leading-power corrections, characterized by a scale larger than \/AMb, can be estimated more reliably, and 
found to be less than 20% of the leading contribution. That is, the heavy-quark expansion makes sense. Note 
that the next-to-leading-power corrections considered here, which can be analyzed under the current knowledge 
of nonperturbative inputs, are not complete. The conclusion drawn in this paper provides a solid theoretical 
base for the PQCD analysis of the A;, baryon charmful decays [43]. 

We have determined the meson wave function from the B — > D^lv decay spectrum, which has a 

maximum at the spectator momentum fraction x ~ 0.36 as expected. This wave function is useful for making 
predictions for the two-body nonlcptonic decays in the PQCD formalism. The results of the B — ► Dir branching 
ratios have been presented in Eq. (69), which are consistent with the experimental data. The detail of this 
subject will be published elsewhere. 
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